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The computer program, MCFIARE, uses Monte Carlo methods to simulate solar 
flare occurrences during an interplanetary space voyage. The total biological 
dose inside a shielded crew compartment due to the flares encountered during 
the voyage is determined. The computer program evaluates the doses obtained 
on a large number of trips having identical trajectories. From these results, 
a dose Dp having a probability p of not being exceeded during the voyage 
can be determined as a function of p for any shield material configuration. 

The user of the code selects any number of solar flares considered to be 
representative of the ones that will occur during future solar active periods 
(these flares are generally selected from the flares that occurred during the 
last solar active period 1956 to 1962). The flares are assumed to occur dur- 
ing these periods. The dose at a distance of 1 AU from the Sun from each of 
these flares behind any shield configuration investigated is input to the 
MCFLARE code. The code accounts for the dependence of the dose received from 
a flare on the distance from the Sun according to a (l/r) a variation, where 
r is the distance from the Sun and the exponent a can be assigned. From 
trajectory parameters, which are input to the computer program, the distance 
from the Sun as a function of time during the trip is calculated. 

To illustrate the use of the code, a trip to Mars and return is calcu- 
lated, and estimated doses behind several thicknesses of aluminum shield and 
water shield are presented. 


INTRODUCTION 

The protons emitted by solar flares are the 
most hazardous source of space radiation encount- 
ered on an interplanetary space voyage which rapid- 
ly traverses the Van Allen belts. To define shield- 
ing requirements for biological protection from 
solar flare protons encountered in manned inter- 
planetary missions, a computer program MCFLARE 
(ref. l) has been prepared which uses Monte Carlo 
methods to simulate solar flare occurrences during 
the trip and records the total associated biological 
doses for the trip behind various shield material 
configurations. 

A set of flare events, considered to be repre- 
sentative of those that will occur during future 
solar active periods, is required. Presently, this 
set is selected from the flares that occurred dur- 
ing the last solar active period from 1956 to 1962. 
The flares are assumed to occur randomly during the 
trip. The doses behind the various shield config- 
urations investigated have to be evaluated for each 
flare, at a distance of 1 AU from the Sun, by using 
another computer code such as one of those described 
in Refs. 2 or 3. (Ref. 4 presents doses from solar 
flare proton spectra behind two shield materials as 
calculated by using Ref. 2. ) These dose values are 
input to the computer program MCFLARE. This pro- 
gram provides for including the dependence of the 
magnitude of the dose received from the flare on 
the distance r from the Sun. This dependence is 
taken to be of the form (r Q /r) a , where r Q is the 
distance between the Sun and the Earth, and the ex- 
ponent a can be any arbitrarily selected value. 

The Monte Carlo sampling procedure (analogous 
to particle transport Monte Carlo) is used to de- 
termine when flares are encountered on the trip, 
which of the input flares occur, and the total dose 
obtained from all flares encountered during the 
trip behind each shield configuration. The comput- 
er program tallies the doses obtained for a large 
number of trips having identical trajectories and 
from this, the distribution of doses incurred on 
the mission for each shield configuration is ob- 


tained. Each dose distribution determines a dose 
D p having the probability p of not being exceed- 
ed as a function of p. 

METHOD OF ANALYSIS 

Any N flare events can be selected as repre- 
sentative of the flare activity that will occur 
during the space voyage. The doses from each flare 
at a distance of 1 AU from the Sun behind one or 
more shield material configurations investigated 
are input to the computer. Each flare is assumed 
to occur randomly on the average of once during the 
solar active period T ref , which is of the order of 
2000 days or 5^ years. 

Monte Carlo Method to Determine Occurrences 
of Flares During Trip 

Inasmuch as flares are assumed to occur ran- 
domly in the solar active period, the probability 
of any flare occurring in the time interval dt is 
(N/T re f)dt = p dt. Then the probability of no 
flare occurring for time t is e - ^ and the 
probability of a flare occurring in dt at time t 
is e^^p dt. 

In the Monte Carlo method, this probability of 
flare occurrence in the time interval dt at time 
t is equated to the probability of selecting a 
random number in an interval about £. If a 

random number set uniformly distributed in the in- 
terval (0,l) is selected, then the probability of 
selecting a random number in d£ about £ is d£. 
Equating these two probabilities, one gets 

e“^p dt = d| 



t =-I ln(l - 4) 
H 
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Inasmuch as both £ 831(1 1 - £ are uniformly — 

distrihuted in the interval (0,l) an equivalent ex- a D = / P^ 1 ~ P I (3) 

pres sion for t which results in the same exponen- P y Mf 2 (D p ) 

tial distribution is 


t = - I In i (1) 

By selecting a random number from a uniformly dis- 
tributed set in the interval (0,l) the time be- 
tween flare events is obtained from Eq. (l). 

When an event has occurred, a particular flare 
of one of the N is also selected by random num- 
bers. The i th flare occurs when a new random 
number occurs in the interval 

< f i = 1, 2, . . N (2) 

N s N 

Dose values input to the code are calculated 
by using the methods of Ref. 2, or 3, for each rep- 
resentative flare using the solar flare proton in- 
tensity and spectral shape as observed at the 
Earth. The code permits an effect of the distance 
from the Sun on the dose (or intensity) received 
from a flare of the form (r 0 /r) a where r Q is the 
distance from the Sun to the Earth, r is the dis- 
tance of the spacecraft from the Sun during the 
flare, and a is an exponent that is assigned. 

The code evaluates r for any time during the 
trip by calculating the mission trajectory from in- 
put parameters describing the two transfer ellipses 
and the total trip time. The first ellipse is the 
trajectory from the Earth to the planet and the 
second describes the return trip. 

The calculation proceeds as follows: Select- 

ing a random number and using Eq. (l) determines 
the time in the mission when a flare occurs. Se- 
lecting another random number and using Eq. (2) de- 
termines the particular flare which occurs. Know- 
ing the time in the mission (and, hence, distance 
from the Sun) and the particular flare, the dose 
received through each shield configuration is 
tallied. The duration of the flare event is then 
added to the time in the mission when the flare 
occurred to give the time when the event is over. 
From a new random number and Eq. (l) , the time 
elapsed until the next flare occurs is calculated, 
and so on. This procedure is repeated until the 
total trip time has elapsed. The total dose from 
all flare events encountered on the trip is deter- 
mined for each of the shields considered. 

Selection of D p 

This procedure is repeated for a large number 
M of trips of identical trajectory. The code di- 
vides the dose range into dose intervals that are 
1 dose unit wide and tabulates the number of trips 
A(D) during which doses that lie within each dose 
interval are encountered. The cumulative fraction 
of the total, trips that encounter doses less than 
the upper dose bound of each interval is also tab- 
ulated (this is an estimate of the probability p 
that a dose equal to the upper dose bound will not 
be exceeded). From this tabulation, for any spec- 
ified p, a dose D p which has the probability p 
of not being exceeded on the trip can be selected 
for any shield. The standard deviation of this 
value from the true value is shown in Ref. 1 to be 


where M is the number of trips in the group in- 
vestigated, and f (D p ) is the dose probability 
density at D p (which is calculated from the tab- 
ulated output) . If A(Dp) is the number of trips 
contained in the interval in which D p occurs, 
f (D p ) is approximately given by 


f(Dp) 


A(D p ) 

M 


If A(D) is a widely fluctuating value in the vi- 
cinity of Dp, averaged values of A(D) can be hand 
calculated from the tabular computer output; 
f(D p ) is obtained from the averaged value of 
A(D p ). Note that for a given p, because f(D p ) 
is independent of the number of trips, 0L p varies 
inversely as >/M. 

To reduce the uncertainty associated with D p 
for a given number of trips, the code is written 
so that the trips constitute a stratified sample 
for the time -to -first -flare selected in accordance 
with the exponential distribution exp[-(Nt/T re f )j . 

ILLUSTRATIVE EXAMPLE 

An interplanetary trip to Mars and return has 
teen evaluated in order to illustrate details of 
the method and to present some representative 
shield requirements for such a mission. Water and 
aluminum were the shield materials considered. 

Selection of Flare Events 

The flare occurrences during the last solar 
active period 1956 to 1962 are assumed to be repre- 
sentative of those that will occur during a future 
active period. Webber in Ref. 5 has compiled a 
record of these flares, their time -integrated 
proton intensities, and spectra. From these rec- 
ords, the 20 largest flares, based on proton inten- 
sity, were selected (the effect of neglecting the 
rest of the flares is small for shields considered 
here). These flares are listed in table I which 
presents the integrated flux of protons having 
energies greater than 30 MeV and some constants A 
and P Q associated with each flare. The spectral 
shape of each flare is assumed to vary as 

. _ . ( Ve 2 + 1876E 1 

N > E = A exp ^ J 

where N > E is the total number of protons per 
square centimeter with energies greater than E. 

The values of P 0 were obtained from Webber's 
compilation, and the values of A were calculated 
to be consistent with flux values of N > 30 MeV. 

It has been observed that some of the flares 
tend to occur in clusters. Table I indicates that 
such clusters occurred in August 1958, July 1959, 
November 1960, and July 1961. These four clusters 
were selected as representative of the clustered 
flare events that may occur and were assumed to be 
events of 12-, 8-, 10-, and 8-day durations, re- 
spectively. The other nine single flares were each 
assumed to have a duration of 2 days. During fu- 
ture solar active periods, each of the four clust- 
ered events and each of the nine single flare 
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events are assumed to occur randomly with a fre- 
quency of once during the active period (taken to 
he 2000 days). 


TABLE I. - FLARES (FROM SOLAR ACTIVE PERIOD 1956 TO 1962) 
SELECTED AS REPRESENTATIVE OF THOSE THAT WILL 


TABLE II. - REM DOSES FROM SELECTED FLARE EVENTS BEHIND 
VARIOUS THICKNESSES OF WATER AND ALUMINUM SHIELD AT 


OCCUR DURING FUTURE ACTIVE PERIODS 


1 ASTRONOMICAL UNIT FROM SUN 


[Integral flux above 30 MeV, and spectral constants A and P for 
each flare are tabulated. ] 


Flare 

date 

N > 30 MeV 
protons/cm 2 

P o 

A, 

protons/cm 2 

2-23-56 

l.O<10 9 

195 

3.4x10® 

1-20-57 

2X10 8 

61 

1. OxlO 9 

8-29-57 

1.2X10 8 

56 

8.6xl0 9 

10-20-57 

5xl0 7 

127 

3.3x10® 

3-23-58 

2.5X10 8 

64 

1. lxlO 9 

7-7-58 

2.5x10® 

62 

1.2X10 9 

8-16-58 

4xl0 7 

64 

1.7X10 9 1 

8-22-58 

7xl0 7 

56 

5. OxlO 9 / Clustered flare event 

8-26-58 

1. lxlO 8 

51 

1.2xl0 10 J 

5-10-59 

9. 6xl0 8 

84 

1.7X10 10 

6-13-59 

8. 5x 10 7 

a 48 

1.2x10*° 

7-10-59 

l.OxlO 9 

104 

1.0X10 10 ! 

7-14-59 

1.3X10 9 

80 

2 6x10 10 f Clustered Dare event 

7-16-59 

9. lxlO 8 

105 

8.9X10 9 J 

9-3-60 

3. 5xl0 7 

127 

2.3X10 8 

11-12-60 

1. 3xl0 9 

124 

8.9xl0 9 | 

11-15-60 

7.2X10 8 

114 

5. 9x 10 9 f Clustered flare event 

11-20-60 

4.5xl0 7 

118 

3.4xl0 9 J 

7-12-61 

4xl0 7 

56 

2 9xl0 9 1 


_ „ fl 


o > Clustered flare event 

7-18-61 

3x10 

102 

3. Ixl0 a J 


Estimated. 


(a) Rem dose behind water shield 


Shield thickness, 

g/cm 2 

Flare date 

1° 

15 

20 

30 

40 

Dose, 

rem 

42. 10 

27.76 

20. 79 

13.03 

9.03 

2-23-56 

1.01 

.38 

. 19 

.07 

.04 

1-20-57 

.43 

.15 

.07 

.03 

.01 

8-29-57 

1.20 

.65 

.42 

.21 

.12 

10-20-57 

1.12 

.40 

.19 

.07 

.04 

3-23-58 

1.26 

.45 

.22 

.08 

.04 

7-7-58 

1.03 

.35 

.17 

.07 

.04 

Aug. 1958 cluster 

11.30 

4.93 

2.67 

1.09 

.57 

5-10-59 

.60 

.20 

.10 

.04 

.02 

6-13-59 

49.69 

23.68 

13. 76 

6.11 

3.36 

July 1959 cluster 

.84 

.46 

.28 

.15 

.09 

9-3-60 

49. 34 

26.49 

16.79 

8.30 

4.85 

Nov. 1960 cluster 

5.48 

2.67 

1. 57 

.71 

.39 

July 1961 cluster 


(b) Rem dose behind aluminum shield 


Shield thickness, 
g/cm 2 

Flare date 

No 

15 

20 

30 

40 

50 

60 

Dose, 

53.93 

37.20 

28. 77 

19.07 

14. 18 

11.15 

9.26 

2.23-56 

2.11 

.86 

.48 

.24 

. 16 

. 12 

. 10 

1-20-57 

.90 

.35 

.19 

.10 

.07 

.06 

.05 

8-29-57 

1.73 

1.00 

.68 

.38 

.25 

.18 

.14 

10-20-57 

2.20 

.89 

.50 

.25 

.17 

.12 

. 11 

3-23-58 

2.48 

1.01 

.56 

.28 

. 19 

. 14 

. 12 

7-7-58 

2.12 

.83 

.45 

.24 

.17 

. 13 

. 11 

Aug. 1958 cluster 

18.98 

9.06 

5. 36 

2.59 

1. 63 

1. 18 

.93 

' 5-10-59 

1.25 

.49 

.26 

.14 

. 10 

.08 

.06 

6-13-59 

78.42 

40. 39 

25. 19 

12.79 

8.12 

5.86 

4. 59 

July 1959 cluster 

1.22 

.71 

1 .48 

.26 

.17 

. 13 

. 10 

9-3-60 

11.89 

41. 17 

; 27. 59 

15.09 

9.87 

7. 17 

5.61 

Nov. 1960 cluster 

8.57 

4.49 

' 2.82 

1. 44 

.92 

.66 

.52 

July 1961 cluster 


Doses For Selected Flare Events 


Discussion of Results 


Table II presents the doses obtained behind 
various shield thicknesses from each of these 13 
events. These doses were calculated by using the 
Lewis Proton Shielding Code, described in Ref. 2, 
and are representative of doses received at the 
center of a spherical crew compartment having the 
given shield thicknesses when the vehicle is lo- 
cated at a distance of 1 AU from the Sun. 

■Table 11(a) presents the rem doses behind the water 
shield, and table Il(b) the rem doses behind the 
aluminum shield. The data presented in this table 
are input to the computer code MCFLARE . 

Trajectory 


Figure 1 shows the distribution of dose ob- 
tained for a large number of identical trips. The 
figure presents the rem dose distribution behind 
20 g per square cm of water shield obtained from a 
computer run of 40 000 trips. The effect of dis- 
tance on dose was assumed to vary as (l/r) 2 . 

In Fig. l(a), A(D), the number of trips on 
which a dose between D - 1 and D has been en- 
countered, is plotted against D. A smoothed 
curve of this distribution A(D) (obtained by av- 
eraging seven values centered about each A(D)) is 
also shown as a dashed line in this figure. 

In Fig. l(b), p(D) the probability of not ex- 
ceeding any dose D is shown plotted against D 


The trajectory selected for the trip is one 
that requires 556 days for the complete mission 
(and happens to be one that results in near mini- 
mum vehicle weight for a departure date in 1983). 
The outward journey from Earth to Mars requires 
280 days, then there is a 40-day stay at Mars, then 
a 236-day return trip during which the vehicle ap- 
proaches within 0.5 AU of the Sun. 


p (D) = i V A(i) 

M f 


i=l 


where M is the number of trips considered in a 
computer run. For the case of Fig. l(b). 


P(D) 


1 

40 000 


D 

£ A(i) 


i=I 
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If Dp is defined as that dose which has a 
probability p of not being exceeded, then from 
Fig. l(b) , Dq . 99 for the group of 40 000 trips, 
occurs in the dose interval between 93 and 94 rem 
and is selected as the upper bound of the interval, 
namely, 94 rem. 

The dose probability density function f(D) is 
approximately equal to A(D)/M. At a Dq. 99 of 
94 rem, the value of f(D), evaluated by using the 
smoothed value of A(D), is 


f(D 0 99 ) = 26 .... = 0.00065/rem 

40 000 


Table III presents values of Do . 99 deter- 
mined by using the computer code for rem doses be- 
hind various thicknesses (expressed as g/cm 2 ) of 
aluminum shield and water shield. These values 
are presented for both the case where the dose re- 
ceived from a flare event is assumed to vary as 
(l/r ) 2 (considered representative of what might be 
a solar proton diffusion model in space) and for 
the case where there is no effect of distance 
(i.e., distance effect, g(r) = l) . Also shown in 
the table are the values of associated 

with each D 0.99 and the corresponding values of 
DO. 99 + 2 °D0.99* 


and the standard deviation from Eq. (3) is 
99 = °* 77 rem * 

There is a 97.7 percent confidence that the 
true dose D g >99 is less than D 0.99 + 2 ODo. 99 ' 
is less than 95.5 rem. 


i j H |i i ! i 


J ssi5i*iP:-5^!s80?'*=3t»===E=.:=- HS s Si§il§H= * *-* - * - = -*:• = 

litilitffillliilllllltHBHillHSilHlHiiHfiHUMiliiMs iliiiHHMHc#; =“ 

■ — ""MiiiuiiiiiiiiUUi! 




i isiji. '.jus; 

UliiSi&alu — - = . . . 

BMIHIHI I IIIIII I II ’ll l U I II I IHir iUJ inmwmnin i n ii mm i inmw^ — 



(A) NUMBER Of TRIPS ENCOUNTER MG DOSES LYWG WfTHM EACH DOSE KTERVAL 


squirt c«ntlmM#r : « 000 trips; dhUnc# MtA wrte as Hi 



TABLE m. - VALUES OF D Q 99 AND CORRESPONDING STANDARD 
DEVIATION BEHIND VARIOUS THICKNESSES OF WATER AND 


ALUMINUM DETERMINED FROM GROUP OF 40 000 TRIPS 


(a) Rem dose behind water shield 


Distance | Quantity 
efl 


Shield thickness, 















±n Fig. 2 the values of Do. 99 + 2 0Do.99 are 
plotted against shield thickness. In Fig.* 2(a), 
the rem dose is plotted against water shield thick- 
ness for cases of g(r) = l/r 2 and g(r) = 1. 
Figure 2(h) is a similar plot of the rem dose be- 
hind the aluminum shield. From curves such as 
Fig. 2 , one can select a shield thickness such that 
on a given mission, there is a probability p of 
not exceeding an accumulated dose Dp. 

These curves indicate how much more effective 
water is as a shield material than aluminum for 
shielding against solar flare protons on the basis 
of grams per square centimeter of material neces- 
sary to maintain a given dose level inside the 
crew compartment. Also shown for this mission is 
that a l/r 2 distance effect on dose received 
from a flare can have an appreciable effect on 
shield requirements when the mission trajectory 
brings the vehicle to within 0.5 AU of the Sun. 

This result indicates a need for information re- 
garding the effect of position from the Sun on ra- 
diation encountered from a flare. 

COMPUTER PROGRAM 

Complete data input instructions for the com- 
puter program MCFLARE are presented in Ref. 1. 

The version of the code presented is operational 
on the Lewis Research Center IBM 7094-11/7044 com- 
puter system. Execution times for the examples 
discussed previously in this report were about 
2 minutes per group of 40 000 trips. The MCFLARE 
code is available from the Radiation Shielding In- 
formation Center of Oak Ridge National Laboratory 
as code package CCC-93. 



<A) REM DOSE BEHIND WATER SHIELD 
Figure 2. - Values of Dq w + 2^ ^ for various 
thicknesses of water and aluminum shields, and 
for cases of distance effect varying as l/r and 
no effect of distance. Number of trips, 40 000. 
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